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O-f. Abstract 

We revisit the old (fourth-order or quadratically generated) gravity model of 
Starobinsky in four space-time dimensions, and derive the (inflaton) scalar po- 
qq ■ tential in the equivalent scalar-tensor gravity model. The inflaton scalar po- 

tential is used to compute the (CMB) observables of inflation, associated with 
curvature perturbations (namely, the scalar and tensor spectral indices, and the 
tensor-to-scalar ratio), including the new next-to-leading- order terms with re- 
spect to the inverse number of e-foldings. The results are compared to the recent 
(WMAP5) experimental bounds. We confirm both mathematical and physical 
equivalence between f(R) gravity theories and the corresponding scalar-tensor 
gravity theories. 
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1 Introduction 



Inflation is a proposal (cosmological paradigm) about the existence of a short 
but fast (exponential, or de-Sitter-type) accelerated grow of the FLRW scale 
factor a(t) in the early Universe, after the Big-Bang but before the radiation- 
dominated epoch [1]. It implies 

a (t) > (1.1) 

Though the whole idea of inflation remains to be a speculation, there is the 
significant (indirect) evidence for it. In the first place, it is the correct prediction 
of CMB fluctuations and large scale structure, in remarkable agreement with 
the WMAP observations of CMB — see eg., ref. [2]. Inflation can generate 
irregularities in the Universe that may lead to the formation of structure. The 
main discriminators among various inflationary models are the spectral indices 
associated with the primordial power spectrum of curvature perturbations [3]. 
For instance, the on-going PLANCK satellite mission is going to provide tight 
constraints on the observable spectral indices with the accuracy of under 0.5 
percent [4]. Though the basic formulae for the spectral indices in terms of 
any inflaton potential are well known [3], their dependence upon the e-foldings 
number can only be computed in a specific inflationary model. Our motivation 
here is to reconsider primary candidates among the inflationary models, as to 
whether they can survive precisional tests in a near future. 

The excellent model of chaotic inflation was proposed by Starobinsky in 
1980 [5]. It is the simplest version of f(R) gravity theories [6], whose extra term 
beyond the standard Einstein-Hilbert term is quadratic in the scalar curvature. 
The Starobinsky model is reviewed in Sec. 2, where we also argue why the 
other (Ricci- and Riemann- curvature) terms in the quadratically generated 
gravitational action are irrelevant to the FLRW dynamics. 

Any f(R) gravity model is known to be mathematically equivalent to the 
certain scalar-tensor gravity via a Legendre-Weyl transform [7] . We review that 
procedure in Sec. 3. However, even in the current literature on the f(R) gravity 
(see ref. [6] and references therein), its physical equivalence to scalar-tensor 
gravity is put into doubt. As is known in Field Theory, any two field theories, 
related by a field redefinition or via duality, have the same observables. In other 
words, the field theories that are mathematically equivalent are also physically 
equivalent. Of course, in specific cases the full equivalence may be very tricky 
(cf., for instance, the AdS/CFT correspondence), so it still makes sense to 
calculate the observables in both equivalent theories. The spectral indices (in 
the leading approximation) of the Starobinsky model were calculated on the 
f(R) gravity side a long time ago [8]. In this paper we do a calculation on the 
corresponding scalar-tensor gravity side. We confirm the leading terms found 
in ref. [8], and calculate the sub- leading corrections to them in Sec. 4, with 
respect to the inverse number of e-foldings. Checking the physical equivalence 
(ie. the same spectral indices) is yet another motivation to our paper. 
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2 Starobinsky model 



There is a priori no reason of restricting the gravitational Lagrangian to the 
standard Einstein-Hilbert term that is linear in the scalar curvature, as long as 
it does not contradict an experiment. The first attempt of that kind was made 
by Weyl as early as 1921. Nowadays, there is no doubt that the extra terms 
of the higher-order in the curvature should appear in the gravitational effective 
action of any Quantum Theory of Gravity. For instance, they do appear in 
String Theory — see eg., ref. [9] for a review. Since the scale of inflation is just 
a few orders less than the Planck scale [3], it is conceivable that the higher- 
order gravitational terms may be instrumental for inflation. It is already the 
case in the simplest modified gravity model having only the terms quadratic in 
the curvature [7]. 

As is well known, there exist only three independent quadratic curvature 
invariants, R^PR^Xp, R^Rfxv and R 2 . In addition, in four space-time di- 
mensions, 

J d A x sf^g (R^R^xp - AR^R^ + R 2 ) (2.2) 
is topological for any metric, whereas 

J d 4 x ^=g~ (m^R^ - R 2 ) (2.3) 

is topological for any FLRW metric. Those combinations do not contribute to 
the (Friedmann) equation of motion for the scale factor, indicating that the 
scalar curvature models play the most important role in cosmological dynam- 
ics. Hence, the most general gravitational action of the highest order 2 in the 
curvature, which may be relevant for inflation, is given by 

S = J d 4 x ^=g~ (2A-R + aR 2 ) (2.4) 

where we have introduced the cosmological constant A and the dimensional 
parameter a = M~ 2 of mass dimension (—2). We use the spacetime signature 
(+, — , — , — ) and the units h = c = 1. The Einstein-Hilbert term in eq. (2.4) has 
the standard normalization with k = Mp^ in terms of the reduced Planck mass 
Mp 2 = 8itGn- The rest of our notation for space-time (Riemann) geometry is 
the same as in ref. [10]. 

The model (2.4) is the simplest representative of the Starobinsky models 
[5]. As was shown in refs. [5, 7], the equations of motion for the action (2.4) 
have an inflationary solution with a/0 (even when A = 0), which is stable 
provided that a > 0. The stability is confirmed by our method in Sec. 3. 

3 f(R) gravity and inflaton 

The model (2.4) is the simplest particular case of the f(R) gravity models 
characterized by an action 

S f = -^j <*V(*) ( 3 - 5 ) 
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with some function f(R) of the scalar curvature. Those models are quite popu- 
lar in the current literature — see eg., the recent reviews [6] and the references 
therein — due to their theoretical applications to inflation and dark energy. 
The gravitational equations of motion derived from the action (3.5) read 

f(R)R^ - \f(R)g^ + g^Uf'(R) - V„V„/'(fl) = (3.6) 

where the primes denote differentiation. Those equations of motion are the 
4th-order differential equations with respect to the metric g^ v (ie. with the 
higher derivatives). Taking the trace of eq. (3.6) yields 

nf(R) + y(R)R-y(R) = o (3.7) 

Hence, in contrast to General Relativity having f'(R) = const., in f(R) gravity 
the field A = f'(R) is dynamical, ie. it represents the independent propagating 
(scalar) degree of freedom. In terms of the fields {g^ u ,A) the equations of 
motion are of the 2nd order in the derivatives of the fields. 

In fact, any f(R) gravity is classically (mathematically) equivalent to a 
scalar-tensor gravity [7] . The equivalence is established by applying a Legendre- 
Weyl transform. The action (3.5) is equivalent to 

S A = j^ J d 4 xV=g {AR - Z(A)} (3.8) 

where the real scalar A(x) is related to the scalar curvature R by the Legendre 
transformation 

R = Z'(A) and f(R) = RA(R) - Z{A(R)) (3.9) 
A Weyl transformation of the metric 

g^{x) (3.10) 
with the arbitrary field parameter 4>{x) yields 

\r - (V^gsTM) « - « 2 <rfyw} 

(3.11) 

Hence, when choosing 

(3.12) 

and ignoring the total derivative, we can rewrite the action (3.8) to the form 

d A x yf^g I ^ + -g^d^dvfy + ^ exp 

(3.13) 

in terms of the physical (and canonically normalized) scalar field 4>{x). 



sv( 3 



exp 



2k4>(x) 
V6 



-gR 



g exp 



2k4>(x) 



A(n(j)) = exp 



-2 K (f)(x) 



4:K(()(X) 
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Equation (3.13) is the standard action of the real dynamical scalar field 4>{x) 
minimally coupled to Einstein gravity and having the scalar potential 



v{4>) 



-2<j) 



.MpiV6. 



(3.14) 



We are now going to employ it as the scalar-tensor gravity model of inflation. 
In order to explicitly derive the inflaton scalar potential (3.14), one has to solve 
for R in terms of (j) by inverting the relation 



f'(R) = A{<t>) 

that follows from eq. (3.9) by differentiation. In the special case of 



f(R) = R-2A- j^R 2 



(3.15) 
(3.16) 



we find 

v(<f>) = 



M ^ 2 + K ) c::p ( 2 ^ } -MhMl ex J^±_\ + MhMl 

8 +A r XP \M Plv / 3j 4 eXP \M Pl V3j 8 

(3.17) 

where the notation A = Mpj A has been introduced. In terms of the new variable 
and the parameter, 



y 



3Mpi 

respectively, the potential (3.17) reads 



and 



V = -M^iM 



v(y) 



v(y) 

V 



A 



l+ -.^-2e* + l 



(3.18) 



(3.19) 



The scalar potential appears to be bounded from below with the only minimum 
at y = (stability!). It is also sufficiently steep for a slow-roll inflation. It is 
the last (third) cosmological term on the right-hand-side of eq. (3.17) that 
dominates in the potential during the slow-roll inflation (when taken alone, it 
gives rise to a de-Sitter inflationary solution), the second term represents the 
lst-order (leading) correction, and the first term is the 2nd-order (subleading) 
correction. 2 In what follows we ignore A. Then the scalar potential for the 
slow-roll inflation gets simplified to 



V(y) = V (e» - 1) 



(3.20) 



A graph of the function v(y) = 
Fig. 1. After a shift <\> — > 4> + 



2y . 



2e y + 1 near its minimum y = is given in 



with 2exp 



2 (f>o 

3 Mpi 



1, the potential (3.20) 



for the sufficiently negative values of y can be approximated as 



V eS (4>) « V 



1 — exp 




(3.21) 



2 As is clear from eq. (3.19), the 'initial' cosmological term A is unimportant during the 
slow-roll inflation. The ratio A/Vb is also negligible from physical (scale) arguments. 
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V 




where we have ignored the subleading contribution. The scalar potential (3.21) 
is known in the inflationary model building [3]. In our treatment of Sec. 4 we 
use the potential (3.20). 

The (R + R 2 ) gravity (or Starobinsky) model is known as the excellent 
model of chaotic inflation in early Universe, and its spectral indices in the 
leading approximation are also known [8]. 3 In the next Sec. 4 we derive those 
indices in the dual (scalar-tensor gravity) picture, and calculate the sub-leading 
terms. 



4 Spectral indices 

The slow-roll inflation parameters are defined by [3] 



and 



V" 
r2 V 



??(0) = Mpj— (4.23) 

where the primes denote the derivatives with respect to the inflaton field 4>. 
A necessary condition for the slow-roll approximation is the smallness of the 
inflation parameters [3], 

e((f>) « 1 and \ v (<j))\ « 1 (4.24) 

The first condition implies eq. (1.1), whereas the second condition guarantees 
that inflation lasts long enough, via domination of the friction term in the 
inflaton equation of motion (in the slow-roll case): 

3H -V' (4.25) 

Here H stands for the Hubble 'constant' H(t) =a /a. Equation (4.25) is to be 
supplemeted by the Friedmann equation 

9 V , 

H = =- (4.26 

3M| V ' 



3 Though it is irrelevant to the early Universe, the Newtonian (weak field) limit of f(R) 
gravity and that of the corresponding scalar-tensor gravity are also the same, as can be easily 
verified by the use of eq. (3.15). 
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It follows from eqs. (4.25) and (4.26) that 



V 

-Mpi-^= < 



(4.27) 



whose solution during the slow-roll inflation (to < ^start < £ < *cnd) is 



^M P1 ln 



_3V3Mpi 



(t - t ) 



(4.28) 



Substituting it into eq. (4.26) and using the definition H =a /a gives rise to a 
differential equation on the scale factor a(t). Its solution is 



a(t) = e Hot 



t-t 
const. 



1-3/4 



(4.29) 



where we have introduced the notation Hp = M/V24. The presence of a 
singularity at t = to in eq. (4.29) is harmless because our inflationary solution 
is only valid during the slow-roll inflation when t > t s tart > to, so that it does 
not apply to the Big Bang. A resolution of the Big Bang singularity is supposed 
to require the higher-order curvature terms in the gravitational effective action 
(2.4). 

The amount of inflation is measured by the e-foldings number 



/•tend i r4> y 



(4.30) 



where the t e nd stands for the (time) end of inflation when one of the slow-roll 
parameters becomes equal to 1. The number of e-foldings between 50 and 100 
is usually considered to be acceptable. 

In the case of the slow-roll inflation with the scalar potential (3.20), we 

find that e{4>) first approaches 1 at cnd = y/fMp^n (2\/3- 3) ps -0.94 M Ph 

since \rj((f))\ approaches 1 later, at </> en( j = — y/|Mpiln| ps —0.62 Mp\. Then 
eq. (4.30) yields 



N e = - {e-y + y) - - (exp 



yf-0.94 - ^| ■ 0.94^ « | (e-» + 3/) - 1.04 



where we have used the notation (3.18). Similarly, we find 

4e 2 y , -4e"(l - 2e^) 



3(l-evy 



and 



7] 



(4.31) 



(4.32) 



Equation (4.31) can now be used to get y in terms of iV e , while a substitution of 
y(N e ) into eq. (4.32) yields both e{N e ) and rj(N e ). The results of our numerical 
calculations (by using MATHEMATICA) are summarized in Table 1. 



7 



An analytic approximation can be obtained by using the expansion with 
respect to the inverse number of e-foldings. For instance, eq. (4.31) yields 



4iV e 16iV e 2 N* V N ! 



9lnN e 0.94 „fln 2 N, 



+ O ~ivr?r~ (4-33) 



Equation (4.32) now implies 



\n 2 N f 



and 

1 31niV e 5 /ln 2 iV e \ fi . 

We are now ready for a calculation of the CMB observable quantitites in 
our inflationary model, ie. for its specific physical predictions. The primordial 
spectrum in the power-law approximation takes the form of k n ~ l in terms of 
the comoving wave number k and the spectral index n. In particular, the slope 
n s of the scalar power spectrum, associated with the density perturbations, is 
given by [3] 

n s = 1 + 2rj - 6e , (4.36) 

the slope of the tensor primordial spectrum, associated with the gravitational 
waves, is given by [3] 

n t = -2s , (4.37) 
whereas the scalar-to-tensor ratio is given by [3] 

r = 16e . (4.38) 

Equations (4.34), (4.35) and (4.36) in our model imply 

2 31niV e 2 /In 2 N e 



The spectral indices are constrained by cosmological observations — see eg. , 
the recent WMAP5 data [11] that implies 

n s = 0.960 ± 0.013 and r < 0.22 (4.40) 

In addition, the amplitude of the initial perturbations, A 2 ^ = M^V / (24ir 2 e) , is 
yet another physical observable, whose experimental value is [3] 

— J = 0.027 M P1 = 6.6 x 10 16 GeV (4.41) 

Equation (4.41) determines the normalization of the i? 2 -term in eq. (2.4) as 
M . /2 _ 2 



4 • J - ■ (2.7) 2 • • 10~ 4 = (3.5 ± 1.2) • 10" 5 (4.42) 



M P1 V 3 v ' ' (i - e y) 2 
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Table 1: The slow-roll parameters and spectral indices for some values of N e 



N e 


e(xl0~ 4 ) 


r)(xW-' 2 ) 


r(xl0~ 3 ) 


n t (xl0- 4 ) 


n s 


35 


5.13 


- 2.56 


8.20 


- 10.3 


0.946 


40 


3.99 


- 2.27 


6.39 


- 7.98 


0.952 


45 


3.20 


- 2.03 


5.12 


- 6.40 


0.957 


50 


2.62 


- 1.84 


4.19 


- 5.24 


0.962 


55 


2.19 


- 1.69 


3.50 


- 4.37 


0.965 


60 


1.85 


- 1.55 


2.96 


- 3.71 


0.968 


65 


1.59 


- 1.44 


2.54 


- 3.18 


0.970 


70 


1.38 


- 1.34 


2.21 


- 2.76 


0.972 


75 


1.21 


- 1.26 


1.93 


- 2.42 


0.974 



where, in the last step, we have used the value of N e = 53.8 ± 18, as it follows 
from eqs. (4.39) and (4.40). The results of our numerical calculations of the 
spectral indices are collected in Table 1. In particular, we find that the WMAP5 
experimental bounds on the scalar spectral index in eq. (4.40) are satisfied in the 
cosmological model (2.4) provided that the e-foldings number N e lies between 
35.9 and 71.8, with the middle value of N e = 53.8. We also find the noticable 
suppression of tensor fluctuations as \r\ < 8.2- 10~ 3 and \nt\ < 10~ 3 . There is a 
possibility of further theoretical modification, which would imply more tuning 
of the spectral indices, when more terms of the higher-order in the curvature 
are added into the action (2.4). 

5 Conclusion 

Our main results are given by eqs. (4.29), (4.31), (4.32), (4.39), (4.42) and Table 
1. The leading terms agree with the known results [8, 12]. We confirm that the 
simplest (Starobinsky) model of (R + R 2 ) gravity with the single new parameter 
M may theoretically describe inflation and still agree with the experimental 
(CMB) observations. As regards the possible extensions to the quartic curvature 
terms, see eg., refs. [13]. The f(R) gravity is extendable to F(R) supergravity 
[14]. 

References 

[1] A.H. Guth, Phys.Rev. D23 (1981) 347; 
A.D. Linde, Phys.Lett. B175 (1986) 395 

[2] D.N. Spergel et al., Astrophys. J. Suppl. 170 (2007) 377 

[3] A.R. Liddle, D.H. Lyth, Cosmological Inflation and Large-scale Structure, 
Cambridge University Press, Cambridge, 2000 

[4] E. Hand, Nature 458 (2009) 820 



9 



[5] A. A. Starobinsky, Phys. Lett. B91 (1980) 99 

[6] S. Capoziello, V. F. Cardone, A. Troisi, Phys. Rev. D71 (2005) 043503, 
astro-ph/0501426; 

S. Nojiri, S.D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4 (2007) 115, 
hep-th/0601213; 

T.P. Sotiriou, V. Faraoni, Rev. Mod. Phys. 82 (2010) 451, arXiv:0805.1726 
[hep-th] ; 

A. De Felice and S. Tsujikawa, f(R) theories, arXiv: 1002.4928 [hep-th] 

[7] B. Whitt, Phys. Lett. B145 (1984) 176; 

V. Miiller, H.-J. Schmidt, A. A. Starobinsky, Phys. Lett. B202 (1988) 198; 
J. D. Barrow, S. Cotsakis, Phys. Lett. B214 (1988) 515; 
K.-I. Maeda, Phys. Rev. D39 (1989) 3159 

[8] V.F. Mukhanov, G.V. Chibisov, JETP Lett. 33 (1981) 532; 
A. A. Starobinsky, Sov. Astron. Lett. 9 (1983) 302 

[9] S.V. Ketov, Quantum Non-linear Sigma-models, Springer- Verlag, Berlin, 
Heidelberg, 2000, Chapters 2,3 and 6 

[10] L.D. Landau, E.M. Lifshitz, The Classical Theory of Fields, Pergamon 
Press, Oxford, 2002 

[11] E. Komatsu et al., Astrophys. J. Suppl. 180 (2009) 330, arXiv:08.03.0547 
[astro-phys] 

[12] A.O. Barvinsky, A.Yu. Kamenshchik, A. A. Starobinsky, JCAP 0811 
(2008) 021, arXiv:0809.2104 [hep-ph] 

[13] M. Iihoshi and S.V. Ketov, Adv. High Energy Phys. (2008) 521389, 
arXiv:0707.3359 [hep-th]; and Int. J. Mod. Phys. A23 (2008) 2153, 
arXiv:0803.2577 [hep-th]; 

S. Kaneda, S.V. Ketov and N. Watanabe, Class, and Quantum Grav. 27 
(2010) 145016, arXiv: 1002.3659 [hep-th]; 

T. Saidov and A. Zhuk, Bouncing inflation in nonlinear (R 2 + R 4 ) gravi- 
tational model, arXiv: 1002.4138 [hep-th] 

[14] S. James Gates, Jr., and S. V. Ketov, Phys. Lett. B674 (2009) 59, 
arXiv:0901.2467[hep-th]; 

S. V. Ketov, Class, and Quantum Grav. 26 (2009) 135006, arXiv:0903.0251 
[hep-th] ; 

S.V. Ketov, F(R) supergravity, arXiv:0910.1165 [hep-th]; in the Proceed- 
ings of the International Conference 'Invisible Universe', Paris, Palais de 
PUNESCO, June 29 - July 3, 2009; http://universe2009.obspm.fr 
S. James Gates, Jr., S. V. Ketov, N. Yunes, Phys. Rev. D80 (2009) 065003, 
arXiv:0906.4978 [hep-th]. 



10 



